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We propose a method to generate a source of spin-polarized cold atoms which are continuously 
extracted and guided from a magneto-optical trap using an atom-diode effect. We show that it 
is possible to create a pipe-like potential by overlapping two optical beams coupled with the two 
transitions of a three-level system in a ladder configuration. With alkali-metal atoms, and in par- 
ticular with ^'^Rb, a proper choice of transitions enables both the potential generation and optical 
pumping, thus polarizing the sample in a given Zeeman state. We extend the Dalibard and Cohen- 
Tannoudji dressed-atom model of radiative forces to the case of a three-level system. We derive 
expressions for the average force and the different sources of momentum diffusion in the resonant, 
non-perturbative regime. We show using numerical simulations that a significant fraction of the 
atoms initially loaded can be guided over several centimeters with output velocities of a few meters 
per second. This would produce a coUimated continuous source of slow spin-polarized atoms suitable 
for atom interferometry. 

PACS numbers: 37.10.Gh, 37.20.-|-j, 32.70.Jz, 05.40.-a 



I. INTRODUCTION 

The achievement of laser cooling [l| resulted in an un- 
precedented control of atomic ensembles, which led to a 
dramatic evolution of atom-based sensors such as atomic 
clocks, inertial sensors, and magnetometers 0- These 
sensors are today state-of-the-art devices and their im- 
provement is a major technological issue. 

However, the cold atom cloud as the source of the sen- 
sors needs to be reloaded after each measurement cycle. 
This limits the maximum repetition rate of the sensor 
and creates a sampling of the measured variable. It is, 
for example, responsible for the Dick effect that is nowa- 
days limiting optical lattice clocks @. The development 
of continuous cold atomic sources can thus lead to impor- 
tant improvements and new designs of atom interferome- 
ters. Continuous atomic sources are often associated with 
an atomic guide that is usually realized with a magnetic 
potential [j, Q and loaded from a cold atomic source. 
Nevertheless, the presence of a magnetic field may be a 
limitation for interferometric applications, and particu- 
larly for magnetometers. 

Here we propose an all-optical method to create a con- 
tinuous source of guided cold atoms loaded directly from 
a magneto-optical trap (MOT), as shown in Fig. [T] (a). 
It relies on an atom-diode effect @, Q that permits the 
atoms to enter into the guide but forbids them to es- 
cape. The method uses optical pumping to both keep the 
atoms in the guided state after loading and spin-polarize 
the sample in a given Zeeman sublevel. The guide thus 
provides a coUimated and continuous source of polarized 
cold atoms suitable for atomic sensors. 

The article analyses the main physical processes in- 
volved in the proposal implementation and is organized 
as follows. In a brief presentation of the method we ex- 
plain how light-shift engineering in a three-level system 
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leads to the continuous extraction of atoms from a MOT 
into an optical guide. We then model the forces for a 
three-level atomic system in a ladder configuration (S- 
system): from the master equation in the Lindblad form, 
we derive the optical Bloch equations and introduce the 
dressed state picture which simplifies the description of 
the internal state dynamics to rate equations. Since the 
guide performance is limited by the heating rates, we 
quantify the various diffusion processes that an atom ex- 
periences when it propagates along the guide. We ex- 
hibit a possible implementation with ®^Rb atoms, but 
the method can be generalized to any alkali-metal atom. 
We finally use the theoretical results previously obtained 
to simulate the stochastic trajectories of atoms in the 
guide solving Langevin-like equations. From these sim- 
ulations we estimate the performance of the guide and 
show that the proposed method is viable. 



II. PRINCIPLE OF THE METHOD 

We now give a brief overview of how the guide is gener- 
ated and continuously loaded, exhibit the various physi- 
cal effects involved, and explain how they combine. 

The main idea is to coaxially overlap two beams 
[Fig. [T] (b)] to engineer the internal state light-shifts and 
design the guiding potential. As shown in Fig. [T] (c) for 
^^Rb atoms, the first beam at 780 nm is tuned to the 
blue of the |5^Si/2,F=2) |5^P3/2) transition, and 

to the red of the |5^Si/2,^' = l) ^ |5^P3/2) transition. 
The 780 nm field then realizes a state-dependent dipole 
potential: repulsive for atoms in \F = 2) and attractive 
for atoms in = 1). The second beam, placed on the 
red of the 5^P3/2 — >■ '^^^5/2 transition at 1529 nm, in- 
duces a differential light-shift on the transition at 780 nm 
0, 1^ . Therefore, the 780 nm radiation is far-detuned at 
the center of the beams, resulting in a strong reduction 
of the dipole potential. As a consequence, the potential 
is created only at the periphery of the beam, produc- 
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FIG. 1. (Color online) (a) Concept: A pipe-like potential is 
generated all-optically. It continuously extracts atoms from 
a magneto-optical trap and guides them, (b) The 780 nm 
and 1529 nm beams are coaxial, with the 780 nm beam be- 
ing slightly larger, (c) Relevant atomic levels of *^Rb ver- 
sus the transverse direction. The 1529 nm beam induces a 
position-dependent differential light-shift on the transition at 
780 nm. The detuning of the 780 nm beam is chosen to gen- 
erate a state-dependent dipole force which implements the 
atom diode. For an atom in \F = 2), the potential in the 
transverse direction is a double barrier, that is, a pipe, due to 
the cylindrical symmetry of the problem. 



ing a state-selective barrier with a pipe shape due to the 
cylindrical symmetry of the problem. Moreover, the dif- 
ferential light-shift also strongly lowers the spontaneous 
emission rate inside the guide, reducing the heating rate 
and thus increasing the guided distance of the atoms. 

Due to the differential light-shift, the 780 nm light 
is resonant (or nearly resonant) with atoms in \F = 1) 
at some position inside the guide. Therefore, atoms in 
\F = 1) from the MOT can pass through the guide barrier 
before being repumped into the barrier-sensitive state 
\F = 2) by the 780 nm light. In other words, an atom 
can enter the guide but cannot leave it: this is the atom- 
diode effect which continuously extracts atoms from the 
MOT and fills the guide. We wiU see in Sec. IVXl that a 
proper choice of polarizations for the optical beams drives 
closed transitions that maintain the atoms in the guided 
state \F = 2) while optically pumping and polarizing the 
atomic sample. 

For the atom-diode effect to occur the 780 nm radiation 
is not far-detuned from the 5^Si/2 ^^^3/2 transition 
and spontaneous emission will play an important role in 
the dynamics of the system. The scattering induces a 



heating of the sample that limits the guided distance for 
a given transport velocity. Moreover, the occupation of 
the excited states is not negligible near the barrier so that 
the dipole force of the 1529 nm light contributes to the 
overall potential. These effects may modify the naive de- 
scription introduced above. Therefore, to verify whether 
the method is viable and what performance to expect, 
we developed a comprehensive model of the problem, 
expanding the Dalibard and Cohen- Tannoudji dressed- 
atom approach . The following sections present a de- 
tailed study of the optical forces and diffusion coefficients 
for a three-level atomic system driven near resonance. 



III. DIPOLE FORCES IN A THREE-LEVEL 
ATOMIC SYSTEM 

We estimate the dipole force and the related momen- 
tum dispersion coefficients in a S-systcm (Fig. [2]). To 
solve this problem, we derive the optical Bloch equations 
for such a configuration and introduce the dressed-atom 
picture, providing a simplified description of the system 
when the secular approximation is valid. 




FIG. 2. Three-level system in a H configuration. The system 
is doubly-driven with electromagnetic fields: the |1) -n- |2) 
transition at ujq'^'' is driven with a detuning Ai = ajp"*^' — uj^j^^ 
(note that with this sign convention A < corresponds to 
blue detuned light) and a Rabi frequency Sli, whereas the 
1 2) o 1 3) transition at cj^^' is driven with a detuning A2 = 
ciJq^'' — uj^j^'' and a Rabi frequency ^2. An atom in |2) can 
decay to |1) at a rate Fi and an atom in |3) can decay to |2) 
at a rate We assume there is no decay from |3) to |1). 



A. Optical Bloch equations 

According to the notations in Fig. [2l the Hamiltonian 
of the atom in the rotating frame is, in the {|1) , |2) , |3)} 
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basis, 

/O \ 

i/at = ^ Ai . (1) 

V A2 y 

By introducing the jump operators ai2 ~ |1) (2| and 
"■23 ~ 1 2) (3|, the atom-hght interaction can be modeled 
with the fohowing Hamiltonian: 

hill ( t \ ^i^2 / t \ 

i/int = ^ ((T12 + Ol^) + ^ + j , (2) 

where the optical phases are assumed to be constant in 
time so that the Rabi frequencies il\ and ^2 are real 
numbers. The Hamiltonian that describes the internal 
energy of an atom coupled to the two optical fields is 
thus V = i/at + i?i„t- 

The internal state evolution of the system is governed 
by the following master equation: 

5tP - \y. p\ + ri£ [ai2] p + V^H [CT23] P, (3) 
Ti 

where the Lindblad superoperators that model the spon- 
taneous emission processes are 

£ [cr] p = crpcr^ - i {crcr^p} , (4) 

and {A, i?} = AB + BA is the anticommutator. Using 
this master equation, the optical Bloch equations for a 
S-system are derived in App. lA II 

In the end, we are interested in the external motion 
of the atom which occurs at timescales much larger than 
the internal state evolution. We thus assume that the 
internal atomic state is always in the steady state poo = 
limi_>oo P- The steady state satisfies dtPoo = and its 
value is derived in App. lA 21 

B. Dressed-atom picture 

The dressed states are defined as the eigcnstates of the 
Hamiltonian V , which admits the following eigendecom- 
position: 

V = U^'VU, (5) 

where V ~ disLg{Ea, Ep, E^) is diagonal in the basis 
{\a) , 1/3) , I7)} called the dressed states basis. U is the 
unitary transformation between the bare and dressed 
states basis {|1) , |2) , |3>} O {|a) , , I7)} M- 

A knowledge of the unitary transformation U allows 
us to define the jump operators for the dressed states as 

= Ua^U^ . Similarly, the density matrix in the dressed 
basis can be obtained from the solution p of the Bloch 
equations in the bare basis: p= UpU^ . 

The main advantage of the dressed-atom description 
is that when the system is strongly driven, the secular 
approximation applies and the equation of motion for the 
internal degrees of freedom reduce to rate equations, as 
shown in App. IB II The decay rate between the dressed 
states |fc) and \k') is 

Tfe^fc, =ri|(5i2),,f + r2|(?23),,f . (6) 



The rate equations are 

dtp + rp = o, (7) 

where p ~ {paa,Pi3i3,Pj'i)'^ is the vector of the dressed 
state populations and the rate matrix is 

r = -r^^,, r^^^+Ff,^^ -r^,^^ . (8) 

The stationary populations p°° are obtained by solving 
the system rp°° = with \\p°°\\^ = 1. 

C. Dipole force 

If the Rabi frequencies ili and are functions of the 
spatial position coordinates r, as expected for Gaussian 
optical beams, then the overall Hamiltonian of the atom, 
including both the internal and external degrees of free- 
dom, is 

H^^ + V{r), (9) 
2m 

where p is the atomic momentum operator and m is the 
mass of a single atom. 

The dipole force operator can be defined from the 
Hcisenberg equation of motion for the momentum op- 
erator [TTI |: 

F = ^P=^[^f,p]=-Vy(r). (10) 

Over a time interval long compared to the lifetimes l/Fi 
and l/r2, the mean force experienced by an atom at 
position r is obtained from the average for the steady 
internal state, 

(F(r)) = Tr [F(r)p^(r)] = -Tr [(Vy(r)) p^(r)] . (11) 

Note that (F(r)) 7^ —V {V{r)) since the density operator 
Poo depends on the position r. 

Since the trace of an operator is independent of the 
basis choice, the mean dipole force can be written using 
the steady state density matrix in the dressed states ba- 
sis: (F) = Tr (Fpoo) = Tr (UFWUpooU^) = Tr(Fpoo) = 
(F). In that case, a dipole force can be associated with 
each dressed state, 

(f) = (f„) + (f^) + (f^), (12) 

where (^F^^ — — p^.VE'fc. The potential can then be 
obtained by integrating the force. 

IV. HEATING RATES 

As we will see in Scc. lVBl the doubly driven S-system 
can create a double barrier potential in exchange for a 
non-negligible spontaneous emission rate near the bar- 
rier position. Those spontaneous emission events are a 
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source of heating. It is thus necessary to quantify the var- 
ious heating processes to estimate whether the proposed 
method is viable. 

The stochastic fluctuations of the atomic momentum 
result in a random walk of the trajectory in momentum 
space that increases the kinetic energy of an atom in the 
guide. This phenomena limits the atomic lifetime and 
thus the guided distance for a given mean velocity of the 
atoms. Two stochastic processes are responsible for this 
effect: 

• the radiation pressure due to the discrete nature of 
the spontaneous scattering events and the random 
direction in which the photon is emitted; 

• the dipole force fluctuation due to the random vari- 
ation of the internal atomic state [III U3 ■ 



The velocity dispersion induced by a given process is 
characterized by a diffusion coefficient D. It is defined 
so that an infinitesimal variation of the velocity during a 
time interval dt satisfies 



(13) 



where dWt = C V dt is the increment of a Wiener noise 
process, and C is a random vector with a normally dis- 
tributed norm of unitary variance {\C\ G A/'(0, 1)) and a 
random direction. 

It must be noted that all the heating processes re- 
sult from the same stochastic source: the spontaneous 
emission. As a consequence, the processes are not only 
correlated but synchronized since the change of the re- 
coil and dipolar forces occurs simultaneously. This effect 
may modify the resulting diffusion coefficient compared 
to the unsynchronized case where the two processes are 
uncorrelated. Nevertheless, in the following we assume 
that the processes are independent and do not consider 
any synchronization effect. 

We will now determine the diffusion coefficients asso- 
ciated with both the radiation pressure and the dipole 
force fluctuation. 



A. Radiation pressure fluctuation 

The diffusion coefficient associated with the radiation 
pressure at the wavelength Xk (fc = 1, 2) is 



(14) 



where kfj is the single photon recoil momentum at A^, m 
is the atomic mass, and the scattering rate is proportional 
to the decay rate and to the population in the bare state 

|fc): rP{r)=r,pr,{r). 

If the diffusion coefScient in a given spatial direction 
is desired, a prefactor must be applied depending on the 
polarization of the considered transition. For example, 
for the isotropic diffusion due to s-wave scattering the 
prefactor is 1/3 for all directions, whereas for a a~^/~ 
transition a factor of 1/4 must be applied for the two 
directions orthogonal to the dipole axis and a factor of 
1/2 for the last direction. 



B. Dipole force fluctuation 

Since the internal atomic state changes at random 
times, the dipole force fluctuates around its mean value 
[Eq. (fT2|) ]. From the rate equations of the dressed states 
populations [Eq. ([7])], we see that the stochastic evolu- 
tion of the dipole force {Ft,t > 0} is a continuous time 
Markovian process with infinitesimal generator T given 
by Eq. dH]). 

The diffusion coefficient is given by Ref. [ll[ : 



t-T t + T 



Ft 



dr. (15) 



Since the process is Markovian, we show in App. lB 2l that 
the diffusion coefficient along a direction e € {x,y,z} 
satisfies the following relation: 



D 



dip 



(16) 



where F^^) = (F^^^ , F^"^ , F^"^ )^ is the vector of the forces 
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in the different dressed states with F^,"^' = — (VE'fe)^, F' is 
the group inverse of F, and the scalar product is defined 
as (x, y) = ^iPu^iUi- This formula allows us to cal- 
culate the diffusion coefficient of any system described 
by a set of rate equations (Markovian process). More- 
over, it provides an efficient way to perform a numerical 
estimation. 



V. IMPLEMENTATION WITH RUBIDIUM 87 
ATOMS 

We now present how to implement the guided contin- 
uous source with ^^Rb atoms. We first show how an ef- 
fective S-system can be realized within the complicated 
structure of ^''Rb transitions using a proper choice of 
polarizations for the optical fields. We then use the pre- 
vious results to determine the parameters for a suitable 
configuration and analyze the atom loading process into 
the guide based on the atom-diode effect. We also ex- 
plain how the transport velocity along the guide can be 
controlled using counter-propagating beams. Finally, we 
simulate atom trajectories along the guide using stochas- 
tic Langevin-like differential equations to quantify the 
expected performance in terms of guided efficiency along 
a given distance and output velocity. 



A. Doubly-closed transitions 

A closed S-system within the ^""Rb structure is real- 
ized using the 5^Si/2 — > 5^P3/2 transition at 780.24 nm 
and the S^Pg/s -> 4^05/2 transition at 1529.37 nm fisj . 
Therefore, it realizes a very good approximation of the 
theoretical model developed above. Note that this idea 
can be generalized to other alkali-metal atoms since they 
have analogous structures. Furthermore, the choice of 
the 1529 nm wavelength is appealing since laser tech- 
nologies developed for optical telecommunications can be 
used. 
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FIG. 3. (Color online) Cycling transitions in **''Rb. (a) The 
use of polarized light allows us to drive a doubly closed 
transition while providing optical pumping, (b) The decay 
from the 4^D5/2 state to the 5^Pi/2 state is forbidden. 



The use of polarized electromagnetic fields 

both optically pumps the sample and drives 
closed transitions between the Zeeman sub- 
states |F 2, rriF = 2) |F' = 3, rn^ = 3) and 
\F' = 3, m'p = 3) ^ \F" = 4, m'^=A), as depicted in 
Fig. [3] (a). Note also that, with this choice of hyper- 
fine states, a decay into the 5^Pi/2 state is forbidden 
{AF > 1), as shown in Fig. [3] (b). Such a configuration 
is thus equivalent to the 5 system sought for. Moreover, 
since an atom in \F ~ 2) loaded inside the guide cannot 
decay into \F ~ 1) after a spontaneous emission event 
the atom will remain in the guided state, thus avoiding 
its loss during propagation. 



B. A possible configuration 

We analyze the properties of the guide obtained in the 
following specific configuration. The 780 nm beam has 
a waist of 1 mm and an intensity of 1000 mW/cm^ (the 
optical power is thus 16 mW), blue detuned by Ai = 
-27r X 485 MHz from the |F = 2) -> |F' = 3) transition 
and Fi = 27r X 6.065 MHz Q. 

The beam at 1529 nm has a waist of 300 /im and thus a 
Rayleigh range of 18 cm. Its intensity is 3 x 10^ mW/cm^ 
(a power of 420 mW) , and the detuning is A2 = +27r x 
133 MHz to the red of the \F' = 3) ^ \F" = 4) tran- 
sition. The transition linewidth is F2 = 27r x 2 MHz 

MM- 

From these parameters, we determine the steady state 
density matrix versus spatial position using the results 
in App. lA 2[ and estimate the dipole potential after inte- 
grating the dipole force [Eq. (fT2)) ]. The results are pre- 
sented in Fig. 131 



C. Loading of the guide 

We see in Fig. |4] (a) that, for an atom in \F = 2), the 
potential has a double barrier shape, creating a pipe-like 
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FIG. 4. (Color online) Dipole trap properties vs transverse 
position at the waist of the beams, (a) Dipole potential for 
an atom in \F = 2) (blue) and \F = 1) (red). The potential 
in IF = 1) is represented only in a region where the scattering 
rate is low since \F — 1) is not a stationary state of the system, 
(b) Scattering rates at 780 nm (blue) and 1529 nm (black) for 
an atom in \F = 2) (solid line) and \F = 1) (dashed line), (c) 
Diffusion coefficient for an atom in \F = 2) due to the recoil 
at 780 nm (blue), at 1529 nm (black) and to the fluctuations 
of the dipole force (brown). 



potential due to the cylindrical symmetry of the prob- 
lem. The guide depth is 365 /iK. Conversely, for an atom 
in \F = 1), the potential is flat at the barrier position. 
The proposed configuration thus implements the state- 
dependent potential necessary to obtain the atom-diode 
effect. 

Nevertheless, for the atom-diode effect to occur an 
atom in \F = 1) must first pass through the barrier and 
then be repumpcd into \F = 2), after passing the barrier 
maximum. In other words, an atom in |_F = 1) should 
not scatter more than a few photons between the mo- 
ment when it enters the 780 nm beam and the moment 
when it reaches the position of the barrier maximum. 
From Fig. [J] (b) , we see that the barrier width is about 
0.1 mm wide and the scattering rate is below 1 kHz before 
the barrier maximum. As a consequence, an atom with a 
velocity of '--^ 10 cm/s scatters less than one photon while 
crossing the barrier, resulting in a low probability for the 
atom to be repumpcd. Conversely, the scattering rate of 
an atom in \F = 1) strongly increases at the guide cen- 
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ter and the atom will be pumped into the trap sensitive 
state, thus realizing the atom-diode based loading. This 
effect is studied in details using a numerical simulation 
in Sec. IVF2l 



D. Heating rates 

By creating a differential light-shift, the 1529 nm radi- 
ation not only creates the guide barriers but also reduces 
the scattering rates, and thus the diffusion coefficients 
due to photon recoil, at 780 nm and 1529 nm by about 
two orders of magnitude [Fig.|3](c)]. This effect strongly 
lowers the heating rate of the guided sample. 

Nevertheless, the heating rate due to the recoils cannot 
be made arbitrarily small. The reason is that the diffu- 
sion due to the dipolc force fluctuation increases with the 
intensity at 1529 nm because of the enhanced light-shift. 
A compromise must thus be found between the spon- 
taneous emission rate and the dipole force diffusion. In 
the present configuration, the parameters are adjusted so 
that the diffusion due to the recoil is of the same order 
as the one induced by the dipole force fluctuation. 



E. Transport along the guide 

To transport the atoms, an extra force must be ap- 
plied along the guide axis. By counterpropagating the 
beams at 780 nm and at 1529 nm, the scattering imbal- 
ance between the two beams pushes the atoms in a given 
direction with a force 

Fp,,h(r) = hdtk = h (r(!)(r)fcW _ (r)4'') u„ 

(17) 

(k) 

where kj^ is the single-photon recoil impulsion at the 
wavelength Afc and is the unit vector in the z direction. 
The spatial distribution of the pushing force is depicted 
in Fig. [51 The velocity of an atom inside the guide is 
almost constant whereas the atom accelerates near the 
barrier. 
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FIG. 5. Pushing acceleration vs the transverse position re- 
sulting from the counterpropagation of the 780 and 1529 nm 
beams {g = 9.81 m/s^). 



F. Numerical simulation 

Since the heating processes limit the lifetime of the 
guided sample, a compromise must be found between 
the transport velocity and the guided length: the atomic 
sample at the guide output should be cold - meaning that 
the output velocity and its dispersion should be of the or- 
der of a few meters per second - while the sample should 
be guided along a few centimeters. This distance is suf- 
flcient, for instance, to deliver atoms inside the magnetic 
shield required for any kind of high-performance atom in- 
terferometer, particularly since the MOT magnetic fleld 
is continuously operated. 

To validate that the proposed conflguration satisflcs 
these requirements we numerically simulate the trajecto- 
ries of atoms along the guide, solving a set of Langevin 
equations. The statistical analysis of the behavior over 
a large number of trajectories allows us to estimate the 
expected performances of the system in terms of guided 
distance and output velocity distribution. Moreover, the 
loading process is included in the simulation to consider 
its influence on the guided sample. 

1. Langevin equations 

The simulation relies on the fact that the internal de- 
grees of freedom of an atom evolve on a much shorter 
time-scale (limited by the lifetime of the excited states to 
submicroseconds) than the external ones (~ ms). This 
allows us to employ a semi-classical description of the 
problem, describing the internal state evolution accord- 
ing to quantum mechanics while modeling the external 
dynamics with Langevin-like equations (Newton equa- 
tions with stochastic velocity increments). 

The momentum diffusion of the atoms generates a ran- 
dom walk of their velocities in the transverse direction, 
inducing stochastic heating. This random walk can be 
modelled from the diffusion coefficients derived in Sec. lIVI 
and with a Wiener increment [Eq. (|13p ] in the differen- 
tial equation of motion. Explicitly, the velocity vector 
V = {vx,Vy,Vz) evolution is governed by the following 
system of stochastic differential equations: 

dv^it) = -F, {r{t)) dt + {r{t))dW^'} 
m 

+ m)dwi'} + ^ D''l^{v{t))dWif\ 

(18a) 

dvy{t) = ^Fy (r(t)) dt + ^D, {v{t))dwl'} 

+ VD2 (r(t))rf<,' + ./Dlilm)dW^T, 

(18b) 

dv,{t) = - [F, {r{t)) + Fpush (r(t))] dt 
m 

+ ^D, (r(t))dW^,(V + VD2 (r(t))dM^i.T 

+ ^/D^m)dW^z'^^\ (18c) 

Wc solve this system using a Monte-Carlo method to 
simulate atom trajectories along the guide and employ 
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FIG. 6. Trajectory of an atom along the guide obtained by numerical integration of the system of stochastic differential 
equations [Eqs. H18|) ]. The simulation includes the loading of the atom inside the guide. The atom is initially in |_F = 1) and 
placed at 1.5 mm from the guide axis. It is launched towards the guide with a velocity distributed according to a Maxwell- 
Boltzmann law with temperature To = 50 /iK. (a) Trajectory along the guide. The atom is reflected off the guide barrier until 
it escapes the guide when the transverse kinetic energy is higher than the height of the barrier, (b) Evolution of the velocities 
along the three directions. When a reflection occurs, the transverse velocities exhibit sudden sign changes and the longitudinal 
velocity increases due to the higher spontaneous emission rate near the barrier, (c) The Brownian motion of the transverse 
velocities responsible for the heating of the sample by increasing the radial kinetic energy. 



the Rungc-Kutta method [17| to perform the numerieal 
integration of the system. 



2. Initial conditions and loading of the guide 

We first describe the initial conditions of the simulation 
and use the simulation to analyze the loading efficiency 

- that is, the probability for an atom to enter the guide 

- and the velocity distribution of the loaded atoms. 
The atom is initially in the \F = 1) state and located 

along the y axis at a distance of 3 mm from the guide axis. 
Since the atoms in the MOT are far from degeneracy, the 
initial velocity vq is distributed according to a Maxwell- 
Boltzmann law 

for a temperature Tq = 50 corresponding to a veloc- 
ity dispersion 5v = yJkBT^^/m ^ 7 cni/s, and the initial 
velocity vector is Vq = (woiO, 0). Moreover, the Doppler 
frequency shift {/S.vd = ^'>^/^) is included in the simula- 
tion. 

During the propagation of the atom towards the guide, 
we estimate at each step of the simulation the proba- 
bility for the atom to be pumped into the stable state 
\F = 2, mp = 4-2) and then randomly determine whether 
the atom is pumped or not. More precisely, we use the 
knowledge of the scattering rate at each position ''(r) 
and the fact that on average 2.4 photons must be scat- 
tered to pump the atom. We consider that an atom is 
loaded inside the guide when it is pumped at less than 
0.5 mm from the guide center (which corresponds to the 



guide radius defined from the position of the barrier max- 
imum). From the simulation of 10"^ stochastic trajecto- 
ries, we find that about 21% of the atoms are loaded into 
the guide. 

After selecting the trajectories that lead to loaded 
atoms, we analyze the velocity distribution of the atoms 
before their further transport along the guide. We define 
the loaded velocity as that when the atom first reaches 
the plane y = 0. From a large number of loading simula- 
tions we estimate the velocity distribution and the result 
is presented in Fig. [71 The resulting distribution is well 
described by the initial Maxwell-Boltzmann distribution 
[Eq. (HH)] for To = 50 /iK shifted by 2.3 cm/s, and the 
mean velocity of the loaded atoms is 13.4 cm/s. This 
deterministic increase of kinetic energy has two contri- 
butions: first the \F = 1) potential is sfightly attractive, 
and secondly the atom rolls down the barrier potential 
after being pumped into \F ~ 2). However, it is inter- 
esting to note that the velocity dispersion of the loaded 
atoms remains that of the atoms in the MOT. 



3. Transport along the guide 

We now analyze the transport of the atoms along the 
guide. An example of a simulated trajectory is presented 
in Fig.[S] We see that the atom is guided along the beam 
axis by reflections off the potential barrier. During the 
reflections, the transverse velocities reverse while the lon- 
gitudinal velocity increases due to the higher scattering 
rate near the guide barrier. This demonstrates the effect 
of the pushing force (discussed in Sec. IV E[) . which is non- 
deterministic since it depends on the stochastic atomic 
trajectory. 
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FIG. 7. (Color online) Velocity distribution of the loaded 
atoms. The histogram is the distribution obtained hom the 
simulation of 2132 Monte-Carlo trajectories, the solid green 
line is the initial Maxwell-Boltzmann distribution with To — 
50 fiK and the solid blue line is the initial distribution shifted 
by 2.3 cm/s. 
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FIG. 8. (Color online) Statistical properties obtained from 
the simulation of 2132 atomic trajectories, (a) Probability 
for an atom to be guided over a given distance. Inset: The 
distribution for the lifetime of an atom in the guide, (b) Lon- 
gitudinal velocity distribution at a guided distance of 5 cm. 



agation distance while reducing the longitudinal velocity, 
or to increase the propagation distance for a given veloc- 
ity. A better control of the longitudinal velocity could 
be achieved by retropropagating the beams at 780 nm 
and 1529 nm and precisely tuning the power ratios be- 
tween the forward and backward directions; in that case, 
a grating would be generated and sub-Doppler cooling 
mechanisms could be exploited. 



VI. CONCLUSION 

We proposed and theoretically analyzed an all-optical 
method to produce a continuous source of cold atoms 
in a spin-polarizing guide. The resulting system is a 
combination of several physical effects, and notably an 
atom-diode effect based on light-shift engineering in a 
three-level atom. 

To engineer the system, we modeled the dipole forces 
in a doubly-driven S-system and provided a deep descrip- 
tion of the diffusion processes responsible for the heating 
of the atomic sample. These results constitute a general 
theoretical framework which can be used to design other 
atom control methods based on S-systems. 

We finally presented and studied in detail a possible 
implementation with ^^Rb. In particular, we isolated 
a closed three-level transition which can be driven with 
reliable optical sources based on diode laser and telecom- 
munication technologies. Since the guide consists in over- 
lapping two optical beams, the experimental implemen- 
tation should be compact and simple to operate which 
are interesting features for embedded applications. We 
showed with numerical simulations that a guiding dis- 
tance of several centimeters is achievable and compatible 
with slow velocities of a few meters per second. 
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Simulating a large number of atomic trajectories allows 
us to estimate the guide performance. The distribution 
of the guided distances - an atom is lost when it is fur- 
ther than 0.5 mm from the guide axis - shows that a 
significant fraction of the loaded atoms (~ 16 %) have 
been guided over at least 5 cm [Fig. [5] (a)]. Moreover, 
from the velocity distribution at the guided distance of 
5 cm [Fig. [5] (b)] we determine that the mean velocity is 
3.9 m/s with a dispersion of 2.1 m/s, corresponding to a 
temperature of about 25 mK. Note that transverse cool- 
ing could be added [Tsl : either to maintain a given prop- 



Appendix A: Optical Bloch equations and steady 
state solution 

1. Optical Bloch equations 

In the {|1) , |2) , |3)} basis, the density operator can be 
written as 

(Pll Pl2 Pl3 \ 
Ph P22 P23 , (Al) 
Ph P*23 P33 J 
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and Eq. ^ provides the optieal Bloch equations for a 
S-system: 



dtPii (pi2 - Pu) + riP22, 

OtP22 = [Pl2 - P12) + (P23 - P23) 

+r2 P33 — riP22, 
19*^33 = (P23 - P23) - r2P33, 

— + lAi I pi2 - i— [P22 - Pll) 



dtPi2 
dtPis. 

dtP23 



+«-^Pl3, 

r2 ^ \ ^2 ^1 

-—+1^2 I Pl3 + *-^Pl2 - «-^P23, 



;(ri+r2)+i(A2-Ai) P23 
(P33 - P22j - «-^Pl3- 



(A2a) 

(A2b) 
(A2c) 

(A2d) 
(A2e) 
(A2f) 



Note that the second equation has been removed since 
it can be obtained from a Unear combination of the first 
and the third equations. 

To solve this system, we spht the density matrix com- 
ponents into real and imaginary parts: pij = Pij + ^Plj- 
From Eqs. (|A3a|) and (jA3b|) we obtain the relationship 
between the coherences and the populations in |2) and 
|3>: 



P22 



P33 



Pl2i 

1 1 

P23- 

t 2 



(A4a) 
(A4b) 



The real parts of Eqs. (jA3c[) . (jA3dp and (jA3e|) provide 
the following relations between the real and imaginary 
parts of the coherence terms: 



2. Steady state of the density matrix 

The steady state solution of the optical Bloch equa- 
tions is obtained for dtp 00 = 0. Since Tr(p) ~ 1, 
one has the following relation between the populations: 
Pll + P22 + y033 = 1- Moreover, pij - p*^ = 2iliapij and 
the steady state thus obeys the following system of equa- 
tions: 



2 



- ililmpi2 + rip22 = 0, 

ri2lm/923 - r2P33 = 0, 

iAi) pi2 - (2p22 + A33 - 1) 



+ «^P13 =0, 

{-^+^/\2)pl3 + ^^Pl2- 



■ni 



P23 = 0, 



r2)+i(A2- Al)] P23 



{P33 ~ P22) - i^-YPi3 = 



(A3a) 
(A3b) 

(A3c) 
(A3d) 

(A3e) 





-^2lTi 
-2A2/r2 



^llT2 



ri+r2 



T1+T2 



P12 

Pl3 
P23 

(A5) 



From this relation and the imaginary parts of Eqs. (jA3c|) . 
()A3dl) . and (|X3e)) . we show that the imaginary parts of 
the coherences are solutions of the following system: 



(A6) 



where 





M 



I Tl+ 4A2 + 2nl + 2^2 (Ai + ^A2 



2no A 



^Ai 



^1^2 1 



4A2 



217i A 



r2 02 



ij^l^l (A2 - Ai) 



(A2 - Ai)j -Of -VLi- 4j^ (A2 - Ai)^ - r2 (Fi + F2) j 



(A7) 



As a consequence, the determination of the inverse ma- 
trix allows us to solve the system (jA6|) . and thus to 
fully determine the steady state density matrix p^o using 
the relations (|X5|) and (|X4)) . 



Appendix B: Calculation of the diffusion coefficient 
of the dipole force fluctuation 



Here we calculate the diffusion coefficient associated 
with the dipole force fluctuation. Using the master equa- 
tion in the dressed state basis we show that, in the secular 
approximation, the dressed state populations obey rate 
equations. This means that the stochastic internal state 
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evolution is a continuous time Markov process. We then 
derive a general expression for the diffusion coefficient 
associated with a Markov process. This result is finally 
applied in the case of the dipole force diffusion. 



1. Evolution of the dressed states populations 

The master equation (U) can be written for the opera- 
tors in the dressed state basis as 



dtp 



V,P 



TiC[a,2]p + r2C[a23]p. (Bl) 



Following [19|, we project Eq. ()Bip on a given dressed 
state Since V \i) = Ei one has {i\ \V ,'p\ \i) — 0, 
the only contribution to the population evolution of the 
dressed states is thus provided by the projection of the 
Lindblad terms: 

{i\ C [a] p\i) = {i\ ap^^ \i) - i {i\ ad^p\i) - i {i\ p^a'' \i) . 

From the decomposition of the density matrix in the 
dressed states basis, 



3,1 



(B3) 



we obtain 



3,1 

We now assume that we are in the limit Vli ^ Fi and 
^2 ^ in which case the coherences evolve fast com- 
pared to the populations. We can thus perform the secu- 
lar approximation where the coherences are replaced by 
their mean value - that is, for j ^ I, (pji) ~ - and thus 
obtain 



I crpo-^ K) 51 Pa f^"^) 



(B5) 



Decomposing again the density operator in the dressed 
state basis, we have 



3,1 



E 



Pj, {i\aa'^ \j 



PU (i| (7(7^ \i) 



(B6) 
(B7) 
(B8) 



where the last line is obtained using the secular approx- 
imation again. From the resolution of unity \ j) (j| = 
1, we finally obtain 



\aa^^i)=p..^\a\Y^\3){]\Y^\i) (B9) 
= P»»Ek'^)yr- (BIO) 



Combining these results according to the master equa- 
tion (jBip , the evolution of the populations finally reduces 
to the following rate equations: 



dtPv, 



where the transition rates are: 

2 



(^12 )i 



F, 



Er. 



(Bll) 



(B12) 



Introducing the dressed states populations vector p = 
{Paa,Pi3l3,Pj'y)'^, thc population equations (|B11|) can be 
written as 



dtp + Tp = 0, 



where 

r = 



3-i-ro 



(B13) 



(B14) 



It is important to note that det F = 0, which means that 
r is singular. 

The populations thus evolve according to: 



p{t) = e-''piO) 



(B15) 



This relation allows us to determine the transition 
probability Pij (r) = P(j,r|i,0), that is, the probability 
to be in the state \j) at time t = t given that it was in 
state \i) at time t = 0. In other words, this is the popu- 
lation Pjjir) given that initially pii{0) = 1. The relation 
(jB15|) thus provides the transition probability matrix 



P(r) 



-,-rT 



(B16) 



The internal state evolution is a homogeneous contin- 
uous time Markov process associated with the transi- 
tion matrix P (r). The Markov process is characterized 
by r which is the infinitesimal generator of this pro- 
cess: it is the generator of the semi-group {P (r) , r > 0} 
and forward time translation is mapped onto this semi- 
group though the Chapman-Kolmogorov relation Vt9 > 
n, P(ri+T2) = P(Ti)P(r2). 



2. Diffusion coefficient of the dipole force 
fluctuation 

We have shown that the stochastic evolution of the 
internal state is a homogeneous continuous time Markov 
process with generator T. Since for each dressed state \k) 
there corresponds a given dipole force Ffe = — Vi?fc, the 
stochastic dipole force evolution {Ft,t > 0} is also a ho- 
mogeneous continuous time Markov process with genera- 
tor r. The transition graph associated with this Markov 
chain is depicted in Fig. [51 

Since there is a unique stationary population distri- 
bution p°°, the Prop. [2] in App. IB 31 gives thc diffusion 
coefficient for the momentum dispersion associated with 
the stochastic evolution of the dipole force. 



^dip — — ^ 



(B17) 



where F — (Fq,F^,F^)^, F" is the group inverse of F, 
and the scalar product is defined as (x, y) = J^iPii^iVi- 
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FIG. 9. (Color online) Transition graph of the Markov chain 
that follows the dipole force in the dressed states basis. 



3. Diffusion coefficient of a Markov process 

Generalizing the definition of the diffusion coefficient 
for the dipole force [Eq. (jl5|) ]. we define here the diffu- 
sion coefficient associated with a homogeneous continu- 
ous time Markov chain (HCTMC) and show that it can 
be efficiently calculated using the group inverse of the 
infinitesimal generator of the process. 

A HCTMC {Xt,t>0} is characterized by its 
transition matrix P (t) with elements (r) = 
P {Xr ~ Xj\XQ ~ Xi). It is a stochastic matrix: 
^jPy (r) = 1, which obeys the Kolmogorov forward 
equation 

-^P (r) + GP (r) = 0, (B18) 

where G is called the infinitesimal generator of the 
HCTMC. Since the initial condition P(0) = 1 must 
be verified, the solution of this differential equation is 
P (r) = exp (— Gr). The generator satisfies J^i Gij = 
and all the off-diagonal elements must be negative. 

The analysis of the long-term behavior of the Markov 
process makes use of the stationary probability vector 
n°° which is the probability to be in a given state of 
the chain after an infinite time. It satisfies Il°°G = 
with ||n°°||i = 1. To quantify the fiuctuations of the 
stochastic variable around the stationary distribution we 
define the diffusion coefficient: 

Definition 1. Let X = {Xt,t>0} be an irreducible 
HCTMC over the countable state space S = {xi} with 
transition matrix P (r) . // a stationary distribution Tl°° 
exists, then the diffusion coefficient of X is defined as 

Dx = ^ [{XooXr) - (Xoo)') dT, (B19) 

where 

{X^Xr) = J2 (r) , (B20) 

ijes 

(Xo,)=^x,nr. (B21) 



The transition matrix elements must satisfy 
limT-_^oo Pij (''■) < oo, which means that any eigen- 
value of G must be positive. Moreover, if all the 
eigenvalues are strictly positive then lim^_j.oo P(''') = 
which is not a stochastic matrix. As a consequence the 
generator must have at least one eigenvalue equal to 
zero, meaning that G is singular. 

These properties - which are a consequence of the 
Perron-Frobenius theorem - allow us to relate the dif- 
fusion coefficient to the generator of the Markov chain. 
More specifically, the diffusion coefficient is the quadratic 
form involving the states of the chain and associated with 
the group inverse (a special kind of pseudoinverse) of G. 
This result is particularly useful for the numerical calcu- 
lation of the diffusion coefficient. 

Proposition 2. Let X ~ {Xt,t > 0} be an irreducible 
HCTMC over the countable state space S ~ {xi} which 
has a stationary distribution 11°° . If G is the infinitesi- 
mal generator of X , then the diffusion coefficient is 

Dx = (x, G«x) , (B22) 

where x = (.t,;) is the column vector of the Markov chain 
states, is the group inverse ofG, and the scalar prod- 
uct is defined as (x, y) = J^ieS^T'^iVi - 

Proof Since {Xoof = E^jgs x^XjUfUf, it follows that 

(B23) 

Moreover, X is irreducible and has a stationary distri- 
bution, and therefore the theorem of convergence to in- 
variant distribution for a continuous time Markov chain 
is satisfied, and Pij{Qo) = lim.r^oo'PijiT) = H°°, which 
results in 

Dx=Yl ^*a;jn°°Zy = (x, Zx), (B24) 

where Z = [Pi''') ^ P(oo)] dr. Finally, Lemma[3]pro- 
vides the relation between the fundamental matrix and 
the generator of the chain: Z = G*. As a result the 
diffusion coefficient is Dx = (x, G'x). □ 

The calculation of the diffusion coefficient results from 
the calculation of the so-called fundamental matrix or 
deviation matrix Z j20l [2]| , which is related to the ex- 
pected time spent in state j starting from i. It is equal 
to the group inverse of the Markov chain generator. 

Lemma 3. 

/>oo 

Zee/ [P(r)-P(oo)] (ir = G«. (B25) 
Jq 



Proof. By the Perron-Frobenius theorem, the generator 
admits the following cigcndecomposition: 

G = V[OeA]V^\ (B26) 
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with A = diag (Ai, . . . , A„), n being the number of 
nonzero eigenvalues and > 0. It resuhs in 



P(t) = e-°" = V [1 e e-^"] V-\ 
and since Ai > 0. 



(B27) 



P(oo) = hm V [1 



-At 



(B28) 



The deviation from the stationary probability is thus 



P(t) - P(oo) = V [O © e-^^] V-\ 
We can now evaluate the integral 

/>oo 



(B29) 



(B30) 



/•OO 

/ V V,, [0 e e-^"] V-j dT (B31) 
Jo ,„ 



-Ar 



dr 



V 



(B32) 



Using again the fact that A.; > 0, it follows that 



e-^^ dr = = diag {X^\ X'^) , (B33) 
and finally 

(B34) 

Writing S = © A and = © A"^ we have 
G = VSV-i and Z = VS»V~i. Clearly SS«S = S, 
S«SS« = and S«S = therefore 

GZG = VSS«SV"i = VSV^i = G, (B35) 
ZGZ = VS«SS'V-i = VS«V-i = Z, (B36) 
GZ = VSS«V-i = VS^SV-i = ZG. (B37) 

In other words, Z has all the properties of the group 
inverse of G, and a matrix that satisfies this properties 
is unique: Z = G". □ 



As an example, we can use this result to calculate the 
diffusion coefficient associated with the dipole force fluc- 
tuation in a two-level atom. 

Example 4. In the case of the dipolar force on a two- 
level system, one as S ^ {+F,~F} where F ~ hVil/2 
and the transitions between the dressed states are de- 
scribed by the following generator fJl /; 



G = r 



cos"' 9 — cos^ 6 
— sin^ 9 sin"' 9 



(B38) 



where the angle 9 satisfies tan 20 = fi/A. The stationary 
population vector is 



n°° = 4fl , (B39) 

where Fpop = F (cos 9 + sin 9j . The generator has the 
following eigendecomposition: G = VAV^^ with A = 
diag(0,rpop) and 



pop 



sm t) cos 
— sin"' 9 sin* i 



The group inverse of the generator is thus 



G» 



V — sin** 9 sin* 

^ pop \ 



(B40) 



(B41) 



The diffusion coefficient is finally obtained from 
Eq. WM) with X = -FY: 



Z?dip = (x, G«x) 



4F2 sin* 9 cos* 9 
r (sin^e + cos"'^)' 



(B42) 



In the high- intensity limit (T2 3> AJ, the angle is 9 = 7r/4 
and the diffusion coefficient becomes -Ddip = fi.^|Vr2p/2r. 
Conversely, in the low-intensity regime (Ti ^ AJ the an- 
gle is small (9^0) and the diffusion coefficient vanishes. 
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